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A note on the fast power series' exponential 

I. S. Sergeev 

Abstract 

It is shown that the exponential of a complex power series up to 



CN 

o 

CN 

C^ \ order n can be implemented via (23/12 + o(l))M(n) binary arithmetic 

operations over C, where M{n) stands for the (smoothed) complexity 
of multiplication of polynomials of degree < n in FFT-model. Yet, 
it is shown how to raise a power series to a constant power with the 
complexity (27/8 + o(l))M(n). 



1 Introduction 

It is very well known that the exponential of a power series (as well as some 
other elementary operations) has the same order of complexity as multi- 
plication, see e.g. [21 [3]. (When we speak about complexity we consider a 
computational model of circuits or straight-line programs over arithmetic 
basis {±, *} U {ax \ a G C}, see e.g. |3j. The field should not be necessary 
complex, it might be real or any algebraically closed, or any other which 

q \ supports appropriate FFT.) 

Previous papers set up a convention to estimate complexity of the basic 
power series' operations (including exponential) in the number of multipli- 
cations of the same size. The very last papers (since 2000) use the special 

/^ '. FFT-way multiplications. 

The complexity function M{n) of FFT-way multiplication can be intro- 
duced as follows. Let F°(n) be the complexity of FFT of order n over C Let 
F{n) be its smoothed version, that is, F(n) = min m >„ F°(m). For simplicity 
we also assume F(n) = Q(n\ogn) (in fact, F(n) = u(n) is sufficient) and 
kF(n) = (1 + o(l))F(kn) for 1 < k < lognQ Then, let M(n) = QF(n). In 
any case M{n) serves as an upper asymptotic estimate of the complexity of 
multiplication of polynomials of degree < n. 



1 These assumptions are for more convenient way of writing final complexity bounds 
only, it does not affect any other aspect of the proof. Notation g = uj(f) means / = o{g). 



Write an upper estimate of the complexity of computing exponent of a 
power series in C[[x]] modulo x n in the form (A + o(l))M(ri). Let us list some 
previously known results: A = 17/6 [lj, A = 14/3 E], A = 13/6 [TO] 
and [I]. Next, we will show that A = 23/12 is also valid. The method is a 
straightforward combination of methods [8j [10] and [I] . To be more precise, 
scheme of computation follows [I] and technique is inherited from [SJ ITU] . 

We will also show that one can compute a constant power of a power series 
in C[[x]] modulo x n with the complexity (B + o(l))M(n), where B = 27/8. 
It slightly improves the previously best known factor B = 41/12 [8j [TO] . 

Some notation. Let / G C[[x]], then /.. n denotes / mod i" and [//^"J 
stands for (/ - f.. n )/x n . If / = J2i>ofi x ^ then A /' J /' ln / ( when /o = x ) 
and e-* (when /o = 0) denotes formal derivative, formal integral, formal 
logarithm and formal exponent respectively: 

t>l i>l i>l i>0 

2 Exponent 

Consider a problem of computing exponent of a power series h, h..\ = 0. 
Denote f = e h ,r = 1/f. Recall that Ah = Af/f. 

The next iterative formula [5] (derived as a solution of an equation Af = 
gf with g = Ah in this case) is valid for m > n: 

f..m+n = f..m + /..„ J (^"V .n [A^..^)/..™/^- 1 ] ) mod X m+ ™. (1) 

Let E(n) and /(n) denote the complexity of computation e^ and 1/f 
modulo x n respectively. Then we can use (1) to compute f m with the com- 
plexity 

E(n) + I{n) + (13 + o(l))F(m) ~ (13/6 + o(l))M (m) (2) 

for appropriately chosen parameters m and n, e.g. n = o(m) and m = 
0(ny/logn). (That is one of the ways to obtain factor 13/6 in the complexity 
estimate for exponent [TO].) 

To achieve the complexity bound (2) split series into blocks of appropriate 
size k, e.g. k e o(n) fl Q(n/y/logn). Then use double DFT of order (2k, k) 
to perform block multiplications in (1). 

Double DFT of order (I, k) as a map from C[x] to C l+k is defined so 
that its first I components are the components of DFT of order I, another k 
components are the components of composition of the variable substitution 



x — > (x and DFT of order k, where ( - - an appropriate complex number. 
Multiple DFT can be defined in similar way [9]. Multiple DFTs are useful to 
perform multiplications of different sizes on the overlapping sets of inputs. 
Double DFT of order (7, k) or its inverse costs as much as DFT of order I 
and DFT of order k plus 0(1 + k) extra operations, see [10] for details. 

In the case I = 2k (as we have) one can use an ordinary DFT of order 3k 
decomposed into outer DFTs of order 3 and inner DFTs of order k instead. 

The main term of the complexity estimate (2) is contributed by: 
3(m/k)F(k) (which we assume to be approximately 3F(m)) operations to 
compute DFTs of blocks of /, the same number of operations to compute 
DFTs of blocks of Ah, the same number of operations to restore blocks of 
the triple product under the integral, 2F(m) operations to compute 2fc-order 
parts of DFTs of blocks of J(. . .) and the same number of operations to re- 
store the product /. n J(. ..). Other steps (precomputation of f,, n and r. n , 
additions, implementation of A and J operators, calculations in the DFT- 
image spaces) contribute o(m log m) in total complexity. 

To provide a hint for verification we consider a subproblem of the triple 
product computation (this step seems to be less evident in the algorithm 
above) in Appendix. Other details (if necessary) see in [10] . 

Next we turn to introduce an improved version of the algorithm with the 
use of idea due to D. Harvey [I]. 

Suppose we are to compute j '„ 2m . Firstly we compute /.. m acting as 
mentioned above. By the way we also have DFTs of blocks of /.. m _ n and 
A{h,, m ) been computed. At the second stage we use formula [2] 

/..2m = f..m + f..m(h - In f.. m )..2m mod X 2 ™ (3) 

derived by the discrete Newton — Raphson method as the solution of an equa- 
tion ip[f] = h with ip[x] = In x in our case. This stage generally consists of 
the two essential parts: computation of In /.. m up to order 2m and final 
multiplication f m by h — lnf m . 

Denote s = A(/ .. m )//.. m . We compute (In/ m ) 2m as Js.. 2m _i via the 
iteration [7] 



.m'+n— 1 S..m'— 1 2" '".. 



Lm/x" 1 '- 1 mod x m ' +n ~\ (4) 



where m' > m. We start from s.. m _i = A(h,, m ). 

To perform the calculations by (4) we use (6+o(l))F(m) extra operations: 
half of them to compute DFTs of remaining blocks of s, another half to 
compute blocks of the triple product (see Appendix for some details). 

To complete the computation of /.2 m we use another (4 + o(l))F(m) 
operations: half of them to compute 2fc-order DFTs of the order of blocks of 



h— In f__ m and roughly the same number to restore blocks of the product in (3). 
Recall that 2fc-order DFTs of almost all blocks of /„ m are also computed at 
the first stage of the algorithm since we use double DFTs (this is the only 
place we gain a benefit from exploiting double DFT). 

To summarize, we can compute the exponent up to order 2ra with the 
complexity (23 + o(l))F(m). 

3 Exponentiation 

Consider a problem of raising of a power series h, h__\ = 1 to a power C € C 
Denote / = h c , r = 1/f, p = 1/h, s = CAh/h. 

The way of computing a power is just similar with that of computing an 
exponent. We will give only a sketch. 

To compute the first half of the required coefficients of / we use the 
formula 

f..m+n = f..m + f..n J (x^V.n [s.. m+n - 1 f.. m /x m ^\ ) mod X m+n (5) 

derived as a solution of the equation Af = sf. Next, we switch to the 
formula 

/..2m = f..m + /..m(Js ~ hi f., m )..2m mod X 2 ™ (6) 

derived as a solution of the equation In / = C In h. 

To solve a subproblem of computing s we use the iteration [7] 

s.. m+n _i = s.. m _i + p.. n {A(h.. m+n ) - s.. m -ih.. m+n ) mod x m+n . (7) 

As above all series are divided into blocks of size k, except that the first 
half of the series Ah is divided into blocks of size 2k. We also use double 
DFTs of order (2k, k). 

Suppose we are given /.. n , r. n , p„ n and s.. n -\. Then we can compute 
/..2m using (40.5 + o(l))F(m) operations. This bound is contributed by the 
following parts: 

(10.5 + o(l))F(m) operations to compute s.. m _i and DFTs of blocks of 
s.. m -i v i a (7). We use DFTs of blocks of h, s, Ah, p and inverse DFTs to 
restore triple multiplications (it is essential that the blocks of A(/i.. m ) are 
double-sized); 

(10 + o(l))F(m) operations to compute another half of s.. 2 m-i- Here each 
iteration (7) is performed via two ordinary multiplications with the use of 
DFTs of order 2k; 

(10+o(l))F(m) operations to compute /. m and DTFs of its blocks by (5). 
The procedure is the same as in the first part of the exponential algorithm; 



(6 + o(l))F(m) operations to compute lia f.. m up to order 2m. This step 
coincides with that of the exponential algorithm above; 

(4 + o(l)).F(m) operations to perform final multiplication in (6) via DFTs 
of order 2k. 

Therefore, we have got the required complexity estimate. 

Research supported in part by RFBR, grants 11-01-00508, 11-01-00792, 
and OMN RAS "Algebraic and combinatorial methods of mathematical cy- 
bernetics and information systems of new generation" program (project 
"Problems of optimal synthesis of control systems" ) . 
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Appendix 

Let f,g,h£ C[[x]]. Consider a problem of computing 

q = f[gh/x m \ mod i", 

where m is a multiple of n. Suppose series f,g,h<E C[[x}} are given divided 
into blocks of size k (we assume for simplicity that n is a multiple of k): 

f = Y1 aixik ' 9 = Y1 hixik > h = Yl Ci%ik ' deg a *' &i ' Ci < k - 

i>0 «>0 i>0 

Suppose we are also given DFTs of order 3k (or, alternatively, double DFTs 
of order (Zi, I2) with l\ +I2 = 3k) of all necessary blocks a^, bi, q. We are to 
show how one can compute q via approximately 3(n/k)F[k) extra operations. 
Flooring makes some complication. We avoid it as following. Let 

Ui= ^ & /^> 9=[u_ 1 /x k \. 

fi+u=m/k+i 



Then 

[gh/x m \ =9 + J2 



mx ik - 



i>0 



Finally we have 



q = 2_] diX tk mod x n , di = ai9 + \. a \ u v 

Note that di are the polynomials of degree < 3k. 

Turn to calculations. Let a* to denote the vector of DFT (or double DFT) 
of polynomial a(x). 

(?) Compute u* for i = —1, . . . , n/k — 1. It costs 0{k)F({m + n)/k) = 
0((m + n) log((m + n) j k)) since u* are the components of the convolution of 
vectors (65, 6J, . . . , &( m+n)/fc _ 1 ) and (eg, c{, . . . , c^^^.J. Recall that 6* and 
c* are the elements of the space C 3k of DFT-images with component-wise 
multiplication. 

(ii) Compute u_i and hence 9 via inverse DFT. It costs 3F(k) + O(k) as 
far as u*_ 1 is known (in fact, 2/c-point inverse DFT suffices here). 

(in) Compute 9*. it costs 3F(k) + O(fc). 



(iv) Compute d* for % = 0,...,n/k — 1. It costs 0(k)F(n/k) = 
0(nlog(n/k)) contributed by a convolution of order n/k in C 3k and 0{n/k) 
extra additions and multiplications in the same space. 

(v) Compute all di and hence q. It costs (n/k)(3F(k) + 0{k)) to compute 
di and 0{n) to restore q. 

Finally we have an upper bound 

3{n/k + 2)F{k) + 0((m + n) log((m + n)/ib)) 
for the total complexity. 



3aMenaHHe o 6bicTpoM BbiHHCJieHHH 
SKcnoHeHTbi CTeneHHoro pa^a 

H. C. CepreeB 

AHHOTail,HH 

rioKa3aHO, hto nepBbie n HjieHOB SKcnoHeHTbi KOMnjieKCHoro CTe- 
neHHoro pa,n,a MoryT 6biTb BbinncjieHbi 3a (23/12 + o(l))M(n) 6imap- 
C"~- ■ Hbix apiicpMeTHHecKHX onepaijHH Hafl C, iyje M(n) o6o3HaHaeT (npn- 

BefleHHyio) cjiojkhoctb yMHOJKeHira MHoronjieHOB CTeneHH < n npn 
noMomn ^1,110. TaKJKe noKa3aHO, hto B03Be^eHne pa^a b nocToaHHyio 
CTeneHb MOJKeT 6bitb peajiH30BaHO co cjiojkhoctbk» (27/8 + o(l))M(n). 

ii- 1 BBe^eHHe 
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Xopoiiio H3BecTHo, hto cjio>KHocTt> BbiHHCJieHHa MHorHx sjieMeHTapHbix one- 



paHHH co CTeneHHbiMH pa^aMH, b tom HHCJie SKcnoHeHTbi, no nopn^Ky Taxaa 



ace, KaK cjiojkhoctb yMHOJKeHHH, cm., HanpHMep, [21 [3]. (Ilofl cjiojkhoctbio 
noHHMaeTcs hhcjio 3jieMeHTOB b cxeMe H3 cpyHKHHOHajibHbix sjieMeHTOB hjih 
HeBeTBHmeitcH nporpaMMe Ha/j; apncpMeTHHecKHM 6a3HcoM {±,*}U{ax | a G 
C}. Bmccto nojia KOMnjieKCHbix HHceji mojkho paccMaTpnBaTb nojie ^chctbh- 
TejibHbix HHceji, npoH3BOJibHoe ajireGpannecKH 3aMKHyToe nojie hjih jno6oe 
jjpyroe nojie, flonycKaiomee ^nO no^xo^amero nopa^xa.) 

IIpe^HiecTByioHiHe paGoTbi cjie^yioT corjianieHHio H3MepaTb cjioacHocrb 
6a30Bbix onepauHft Ha,n CTeneHHbiMH pajj;aMH (BKjnonaH 3KcnoHeHTy) hhcjiom 
yMHoaceHHH 3KBHBajieHTHoro pa3Mepa. B nocjie^HHx paGoTax (nocjie 2000 r.) 
Hcnojib3yioTCH cneuHajibHbie yMHoaceHHH, ocHOBaHHbie Ha ^HO. 

OyHKHHK) CJIOJKHOCTH M(n) ffil O-MeTO^a yMHOJKeHHH MOJKHO BBeCTH 

cjie/jyioiHHM o6pa30M. 06o3HanHM nepe3 F°(n) cjiojkhoctb ,3,110 nopaflica 
n Hajo; C, a nepe3 F(n) — npHBejo;eHHyio cjiojkhoctb, F(n) = min m > n F°(m). 
npe^nojioacHM jj;jia yn;o6cTBa, hto F(n) = Q(nlogn) (Ha caMOM jj;ejie, jj;o- 
CTaTOHHo 6bijio 6m ycjioBHH F(n) = w(n)) h kF(n) = (1 + o(l))F(kn) npn 
1 < k < logn. (yKa3aHHbie npe^nojiojKeHHH Hcnojib3yiOTCH tojibko jj;jih 
yn,o6cTBa 3anncH nojiyneHHbix oh^hok cjiojkhocth h He 3aTparHBaioT ^pyrne 
acneKTbi jj;oKa3aTejibCTBa. 3anncb g = oo(f) o3HanaeT / = o(g).) IIojiojkhm, 



M(n) = 6F(n). B jiio6om cjiynae M(n) cjiyacnT acnMnTOTnnecKon BepxHeft 
on,eHKon cjiojkhocth yMHoaceHna MHoronjieHOB CTeneHH < n. 

3anHineM Bepxmoio ouemcy cjiojkhocth BBinncjieHna: SKcnoHeHTBi kom- 
njieKCHoro cTeneHHoro pa^a no Mo^yjiio x n b BH^e (A + o(l))M(n) h ne- 
penncjinM H3BecTHbie paHee pe3yjibTaTbi: A = 17/6 [I], A = 14/3 [3 [6], 
A = 13/6 [HI |TU] n |3]. HnjKe 6yn;eT noKa3aHo, hto on,eHKa cnpaBe^jiHBa npn 
A = 23/12. IIpefljiaraeMbiH MeTo,n; nojiynaeTca Henocpe^cTBeHHtiM kom6h- 
HHpoBaHHeM MeTo^oB |S1 [10] n pQ. Bojiee tohho, cxeMa BbinncjieHnn cjie- 
pyeT ajiropHTMy jJJ, ho npn stom ncnojiB3yiOTCs TexHnnecKne npneMbi H3 
pa6oT [El [TDJ. 

TaKace 6yn;eT noxasaHo, nro B03Be^eHHe KOMnjiexcHoro CTeneHHoro pa^a 
b nocTOHHHyio cTeneHb no Mo^yjno x n mojkho peajiH30BaTb co cjiojkhoctbio 
(B + o(l))M(n), rflfi B = 27/8. 3to nyTb jiynnie paHee nojiyneHHofl oneHKn 
B = 41/12 [HI [10]. 

BBe^eM Heo6xo^HMBie o6o3HaneHna. LlycTb / e C [[#]]. Tor^a nepe3 f__ n 
6yn;eM o6o3HanaTB / mod x n , a nepe3 \_f/x n \ cooTBeTCTBeHHo (/ — f„ n )/x n . 

Ecjih / = J2i>ofi x \ TO A /> J /> ln / ( n P H /o = 1) h e-/ 1 (npn / = 0) 
o6o3HanaioT (popMajiBHbie npoH3BO^Hyio, HHTerpaji, jiorapncpM n SKcnoHeHTy 

COOTBeTCTBeHHO : 

i>l i>l i>l i>0 

2 9KcnoHeHTa 

PaccMOTpnM 3a^any BBinncjieHna 3KcnoHeHTBi CTeneHHoro pa^a h, h„i = 0. 
06o3HannM / = e h , r = 1/ f. 3aMeTHM, hto Ah = Af/f. 

Cjie^yioinaH HTepannoHHaa: cpopMyjia [5] (nojiynaeMaa H3 ypaBHeHna bh- 
,n,a Af = gf, rflfi b HanieM cjiynae g = Ah) cnpaBe^jiHBa npn m > n: 

f..m+n = f..m + /..„ J {x^Y.n [A(h . , m+n ) f . , m / X^ 1 \ ) mod X m+n . (1) 

06o3HannM nepe3 E{n) n I(n) cjiojkhoctb BbinncjieHna e^ n 1/f no mo- 
flyjno a; n cooTBeTCTBeHHo. Tor^a, ncnojib3ya (1), mojkho BBinncjinTb f\ m co 

CJIOJKHOCTBK) 

E(n) + I{n) + (13 + o(l))F(m) ~ (13/6 + o(l))M(m), (2) 

ecjin napaMeTpbi m n n Bbi6paTb noflxoflamHM o6pa30M, HanpnMep, n = 
o(m) n m = 0(n\/Iogn). (3to oahh H3 cnocoGoB bbihhcjihtb 3KcnoHeHTy co 

CJIOJKHOCTbK) 13/6 yMHOJKeHHH [10 J.) 



Jliis ^0Ka3aTejibCTBa oijeHKH (2) pa3o6teM paflM Ha Gjiokh no/rxo/XHinero 
pa3Mepa k, HanpnMep, k G o(n) fl Q(n/ y/log ri). JIjik BBinojmeHns Gjiohhbix 
yMHoaceHHit b (1) 6yn;eM Hcnojib30BaTb /xBoftHBie ,11,110 nopa/xita (2k, k). 

fleounoe JJI7<P nopa^Ka (I, k) kok oToGpaaceHne H3 C[x] b C l+k onpe^ejia- 
eTca cjie/xyioiirHM o6pa30M: nepBbie I KOMnoHeHT hbjijiiotch KOMnoHeHTaMn 
^1,110 nopa^Ka /, jxpyrne k KOMnoHeHT nojiynaiOTca: KOMno3HHiieH 3aMeHBi 
nepeMeHHofl x — > (x h ,1111 <J> nopa/rica k, r^e ( -- no/xxoiXiiinee KOMnjieKC- 
Hoe hhcjio. IIo^o6hbim o6pa30M onpe^ejiaioTca ,2,11(1? GojiBinefl KpaTHocTH. 
MHoroxpaTHBie ,U,n<J> BBiro/XHo HcnojiB30BaTb b cjiynae, Kor^a BbinojmaeTca 

HeCKOJIbKO yMHOJKeHHfl pa3HOfl pa3MepHOCTH, HO C oGlHHMH BXO^aMH. ,HbOH- 

Hoe ,H,n<J> nopa^xa (l,k), paBHo xax h o6paTHoe k HeMy, peajiH3yeTca co 
cjioJKHOCTbio /^nO nopa/iKa I, ,2,11 <J> nopa/XKa k njnoc 0(1 + k) ^onojiHH- 
TejibHbix onepanrni, no/xpoGHee cm. b [TU] (o MHoroxpaTHbix ,2,11 <J> cm. [9j). 

B cjiynae / = 2k (T.e. b HaineM cjiynae) BMecTo ^bohhbix ,2110 mojkho 
HcnojiB30BaTb o6biHHbie ^1,110 nopa/XKa 3k, npeflCTaBjiaeMBie KOMno3HiriieH 
BHeniHHx ^1,110 nopn^Ka 3 h BHyTpeHHHx 2420 nopa/jKa k. 

rjiaBHbiH HjieH b oueHKe (2) cKjia/XBiBaeTca: H3 cjie/xyioiiiHx cocTaBjiaio- 
uihx: 3(m/k)F(k) (hto no npeflnojKxaceHHio npH6jiH3HTejibHO paBHo 3F(m)) 
onepau;HH rjik BBinncjieHHH 2420 Gjiokob pa^a /, ctojibko ace — fljra bbi- 

HHCJieHHH ^nO 6JIOKOB pa^a Ah, CTOJIBKO >Ke — flJIH BOCCTaHOBJieHHH 6jio- 

kob TpofiHoro npoH3Bejo;eHH2 nojj; 3HaKOM HHTerpajia, 2F(m) onepaniiH fljia 
BBiHHCJieHHH 242 ^ nopa/iKa 2k Gjiokob psija;a J(. . .) h ctojibko >Ke — jjjtsi pe- 
KOHCTpyKiniH npoH3Beji;eHH2 /.. n J(. ..). OerajiBHBie ji;eHCTBHH (npejxBBinnc- 
jieHHe /.. n h r.. n , peajiH3anH5i onepaTopoB A h J, BBiHHCJiefflra b npocTpaHCTBe 
OypBe-o6pa30B) bhocst bkji&ji o(m log m) b o6myio oueHKy cjiojkhocth. 

C n;ejiBio npeflocTaBHTB HeKOTopoe o6ocHOBaHHe yKa3aHHBiM oireHKaM b 
IIpHjiojKeHHH nojxpo6Hee paccMaTpHBaeTca no/x^a/rana TpoflHoro yMHoaceHHa 
(BepoaTHO, HaHMeHee oneBH^Haa nacTB onncaHHoro ajiropHTMa). 2,pyrne Jie- 
TajiH ji;oKa3aTejiBCTBa MoryT 6bitb BoccTaHOBjieHBi no paGoTe [IIP] . 

TenepB h3Jio>khm yjiynnieHHyio Bepcnio ajiropHTMa, ncnojiB3y5i n^eio H3 
pa6oTBi 2,. XapBH [4] . 

IlycTB Tpe6yeTca bbihhcjihtb f,,2m- A 713 SToro CHanajia mbi BbiHHCJiaeM 
/.. m cnocoGoM, yKa3aHHBiM BBinie. Ilpn stom nonyTHO TaKace bbihhcjihiotcsi 
2420 6jiokob f\, m - n H A(h„ m ). Ha btopom 3Tane ncnojiB3yeTCH cpopMyjia [2] 

f..2m = f..m + f..m(h-\nf„ m )„2m mod X 2 ™ , (3) 

nojiynaeMaa jxhckpcthbim aHajioroM MeTOJi,a HbioTOHa — PacpcoHa H3 ypaB- 
HeHHii ip[f] = h, rp<5 b HanieM cjiynae ip[x] = lnx. 3Tan coctoht H3 ppyx 
ochobhbix nacTeii: BBinncjieHne In /.. m no MOji;yjiio x 2m n 3aKjnonnTejiBHoe 
yMHoaceHHe /.. m Ha h — In f\ m . 



06o3HaHHM s = A(/ ..m)/ f..m- Tor,na (ln/.. m )..2 m BbmncjiaeTca kbk Js..2m-i 
npn noMoniH HTepaHHOHHoii cpopMyjibi [7J 



.m'+n— 1 °..m'— 1 •*' ^.. 



m'-i/.m/x™^ modx^"- 1 , (4) 



r,ne m' > m. Hanajio pa^a s H3BecTHo: s.. m _i = A(/i.. m ). 

JXjis BbiHHCJieHHH no cpopMyjie (4) Hcnojib3yeTCH jj;onojiHHTejibHo (6 + 
o(l)) F(m) onepannft: nojiOBHHa H3 hhx — pjm BbiHHCJieHHH HefloCTaromHX 
,11,110 Gjiokob pa^a s, jjpyran nojiOBHHa — jj;jih BoccTaHOBjieHHH Gjiokob Tpoii- 

HOrO npOH3Be^eHHH (o TpOHHOM yMHOJKeHHH CM. B IIpHJIOJKeHHH) . 

JJjis. 3aBepineHHa BbiHHCJieHHH /..2m BbinojmHeTCH eme (4 + o(l))F(jn) 
onepanHft: nojiOBHHa H3 hhx — jj;jih BbiHHCjieHHH ^nO nopHijKa 2k 6jiokob 
pa^a h—lnf__ m , npHMepHo cTOJibKo >Ke — fljis BoccTaHOBjieHHH 6jiokob npoH3- 
Be^eHHa b cpopMyjie (3). 3aMeTHM, hto ,3,11 <J> nopH/jKa 2k noHTH Bcex Gjiokob 
f__ m BBiHHCJieHti Ha nepBOM 3Tane ajiropHTMa xax nacTH jtbokhbix ,3,11(1? (3to 
to MecTO, r^e npHMeHeHHe ^bohhbix ,3,11 <J> ^aeT Bbmrpbnn). 

OKOHHaTejibHo fljiH BBiHHCJieHHH nepBbix 2m hjichob pa^a 3KcnoHeHTbi 
Mbi nojiynaeM ou,eHKy cjiojkhocth (23 + o(l))F(m). 

3 BosBe^eHne b CTeneHb 

PaccMOTpHM 3a^aHy B03Be/i;eHHH pa^a h, h..i = 1 b CTeneHb C G C 06o3Ha- 
™ f = h c ,r = 1/f, p = 1/h, s = CAh/h. 

JXjis B03Be^eHHH b CTeneHb mojkho npe^Jio>KHTb ajiropHTM, aHajiorHH- 
Hbifi onncaHHOMy Bbirne ajiropHTMy BbiHHCJieHHH SKcnoHeHTbi. H3Jio>khm ero 

KpaTKO. 

Mjia^niaa nojiOBHHa hckombix KoscpHinieHTOB pafla / BbiHHCJiHeTCH no 
cpopMyjie 

J..m+n = J..m + J..n J [% r ..n [_ s ..m+n-lj ..m/ x J) mod X , (5) 

ocHOBaHHon Ha cooTHonieHHH A/ = sf. ^ajibHenniHe BbiHHCjieHHH Bbinoji- 
hhiotch no cpopMyjie 

/..2m = /.m + f..m(Js - In /. m )..2m mod X 2 ™ , (6) 

KOTopan bbibojjhtch H3 ypaBHeHHH In / = C In h. 

fljisi BbiHHCJieHHH pima s Hcnojib3yeTC2 HTepan;HOHHaH cpopMyrra [7J 

s ..m+n-l = s ..m-l + P..n(A(/l.. m +n) — s ..m-l' l ..m+n) HlOCl X . (7) 



Kax h npeac^e, Bee pa^bi pasGnBaioTca Ha Gjiokh ^jihhbi k, 3a hckjho- 
neHHeM nepBofi nojioBHHbi pa^a Ah, KOTopaa pa36HBaeTcs Ha Gjiokh ajikhbi 
2k. 3aTeM b BbiuncjieHHax Hcnojib3yioTca: ^BoitHbie ,11,110 nopaflica (2k, k). 

IlycTb f__ n , r„ n , p.. n h s.. n _i flaHH. Tor^a mbi MoaceM bbihhcjihtb /..2m sa 
(40,5 + o(l))F(rn) ^onojiHHTejiBHBix onepanHii. 3Ta orjemca cKjia^tiBaeTca 
H3 cjie/ryiomHX cjiaraeMbix: 

(10,5 + o(l))F(m) onepaiiHH fljia: BbiHHcneHHa s.. m _i h ,11110 6jiokob 
s„m-i n P H noMoriiH (7). A- 713 btoto bbihhcji3iotc2 ,11110 6jiokob ph^ob h, s, 
Ah, p h o6paTHbie ,U,nO rjik peKOHCTpyKHHH Tpofmbix npoH3Be^eHHH (oueH- 
Ka ^ocTHraeTca Gjiaro^apa TOMy, hto Gjiokh psma A(/i.. m ) hmciot ^bohhoh 
pa3Mep); 

(10 + o(l))F(m) onepaHHH fljis BbiuHcneHHS CTapniefl nojioBHHbi s..2 TO -i- 
3^ecb Kaac^aa HTepairna: no cpopMyjie (7) BbinojmaeTCH 3a ^Ba oGbiHHbix 
yMHoaceHHa, jj;jia nero Hcnojib3yiOTca: ,11,110 nopa^Ka 2k; 

(10+o(l))F(m) onepauHH ^jia BbiHHCjieHHa uacruuHoro pa^a /.. m h ,3,140 
ero 6jiokob (5). I4poH,e,nypa Taxaa ace, Kax b nepBoft Hacru ajiropHTMa bw- 
HHCJieHHa SKcnoHeHTbi; 

(6 + o(l))F(m) onepauHH fljia BbiHHCJieHHH 2m nepBbix hjichob pa^a 
ln/.. m . 3tot niar coBna^aeT c cooTBeTCTByioiiiHM niaroM onncaHHoro Bbinie 
ajiropHTMa fljis scnoHeHTbi; 

(4+o(l))F(m) onepau,HH jj;jia peajiH3aiiHH 3aKjnouHTejibHoro yMHoxceHHa 
no cpopMyjie (6) npn noMomH ,H,I40 nopa^Ka 2k. 

CyMMnpya, nojiynaeM 3aHBjieHHyio oueHKy. 

Pa6oTa BbinojiHeHa npn (pHHaHCOBoit no^mepacxe POOH, npoeKTbi 11- 
01-00508 n 11-01-00792-a, n nporpaMMbi cpyH^aMeHTajibHbix Hccjie^oBa- 
hhh OMH PAH «Ajire6pannecKne n KOM6nHaTopHbie MeTo^bi MaTeMaTHue- 
CKoit Kn6epHeTHKH n HHcbopMaHHOHHbie cncTeMbi HOBoro noKOJieHHH» (npo- 
eKT «3a^ann onTHMajibHoro cnnresa ynpaBjiaioninx cncTeM»). 
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IlpHjic>;>KeHHe 

nycTB f,g,h£ C[[x\], h Tpe6yeTca bbihhcjihtb 

q = f[gh/x m \ mod i", 

r,n;e n\m. nycTB pa^Bi f,g,h G C[[x]] pa36nTBi Ha Gjiokh ajihhbi k (fljia 
npocTOTbi Mbi nojiaraeM, hto n ^cjihtch Ha k): 

f = Yl ai%ik ' 9 = Y1 hi%ik i h = Yl Ci%ik ' deg a " bi > Ci < k - 

i>0 i>0 i>0 

nycTB TaKace ^aHBi Beicropa /Hi® nopa^Ka 3k (hjih ^bohhbix /Hi® nopa^xa 
{hjh), r^;e l\ + l 2 = 3k) Bcex Heo6xo,a;HMi>ix 6jiokob Oj, bi, Ci. noKaaceM, 

HTO q MOJKHO BBIHHCJIHTB npH6jIH3HTejIBHO 3a 3(n/k)F(k) ^OnOJIHHTejIBHBIX 

onepaHHii. 



HeKOTopbie TpynHoc™ co3^aeT oTGpacbmaHHe Mjia^innx hjichob pa^a gh. 
Hx mojkho npequojieTb cjie^yiomHM o6pa30M. 06o3HaHHM 

Ui= ^ h v c v> 0=[u_ 1 /x k \. 

fi+v=m/k+i 



Tor^a 



[gh/x m \ =9 + J2^ xlk - 



i>0 



OKOHHaTejibHo nojiynaeM 



q = 2_. diX tk mod x n , di — ai9 + \. a \ u ^- 

3aMeTHM, hto di — MHoronjieHbi cTeneHH < 3k. 

OnnineM nopa^oK BbiHucjiemift. Hepe3 a* 6yn;eM oGosHanaTb BeKTop ,2,11 <J> 
(hjih ^BOHHoro ,H,nO) MHoroHjieHa a(x). 

(i) Bbihhcjihiotch u*, r^e % = — 1, . . . ,n/k — 1. Cjiojkhoctb mojkho on;e- 
hhtb Kax 0(k)F((m + n)/k) = 0((m + n) log((m + n)/k)), t.k. u* — kom- 

HOHeHTBI CBepTKH BeKTOpOB (6*, 6*, . . . , 6( m+n)/fc _ 1 ) H (eg, c\, . . . , C^)/^). 

HanoMHHM, hto b* h c* npHHa^jiexcaT npocTpaHCTBy <J>ypbe-o6pa30B C 3fe c 
noKOMnoHeHTHofi onepauHeii yMHcoKemra. 

(ii) BBiHHCJiaeTca u_\ h, cjie^oBaTejiBHo, # npn noMouiH o6paTHoro ^nO. 
Cjio>KHOCTb oijemiBaeTCJi Kax 3F(k) + O(k), t.k. BeKTop u!_ 1 H3BecTeH (b 
^eitcTBHTejibHocTH, ^ocTaTOHHo BbinojiHHTb oGpaTHoe An® nop$mKa 2k). 

(Hi) BbiHHCJiHeTCH $*. Cjio>KHocTb 3F(k) + O(k). 

(w) BblHHCJIHHJTCH d* , Fflfi 1 = 0,..., fl/k — 1. CjICOKHOCTb 0(k)F(n/k) = 

0(nlog(n/k)) cKjiaflbraaeTCH H3 cjicokhocth CBepTKH nopa^Ka n/k b C 3fc h 
0(n/k) ^onojiHHTejibHbix cjioaceHufi: b tom :»ce npocTpaHCTBe. 

(v) BbmncjiHioTCH Bee d{ h 3aTeM g. 3,n;ecb (n//c)(3F(A;) + 0(/c)) onepamrft 
Hcnojib3yeTca ,zyiH BbiHucjiemiH di h O(n) — fljra BoccTaHOBjieHiis g. 

OKOHHaTejibHo ,zyiH ajiropHTMa b u,ejioM nojiynaeM orjeHKy cjicokhocth 

3{n/k + 2)F(k) + 0((m + n) log((m + n)/ife)). 



